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FEASIBILITY OF DETERMINING ATMOSPHERIC OZONE 
FROM OUTGOING INFRARED ENERGY 

C. Prabhakara 
ABSTRACT 

The high resolution (5 cm -1 ) measurements of the outgoing in- 
frared energy in the region of the 9.6 p. ozone band offer a means of 
determining the vertical distribution and total amount of ozone in 
the earth's atmosphere. With the application of radiative transfer 
theory and perturbation technique a method is developed to deduce 
such information. The method hinges on a two parametric repre- 
sentation of the ozone distribution in the earth's atmosphere. 

An error analysis based on four case studies is presented to 
show how well the atmospheric ozone could be determined. It is 
found that a small error in radiance value is magnified consid- 
erably in the inferred atmospheric ozone. 
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INTRODUCTION 

Atmospheric ozone is photochemieally produced by solar near ultraviolet 
radiation. The global distribution of this gas as a function of latitude and season 
can reveal the stratospheric circulation, and coupling between the troposphere 
and stratosphere. Day to day changes in the weather are also reflected in the 
ozone variations. Thus it appears that a method of indirectly obtaining the at- 
mospheric ozone from satellite borne radiation sensors could be very useful. 

The problem of indirectly determining the vertical distribution of tempera- 
ture and water vapor in the earth's atmosphere from remote radiation sensors 
has been studied by several investigators (see Wark and Fleming, 1966 and also 
for a review see Conrath, 1968). Such investigations were concerned with the 
infrared spectral regions of the 15/a C0 2 band, and the water vapor 6.3/a vibration- 
rotation and 20p. rotation bands. However in addition to these bands the IRIS 
(Infrared Interferometer Spectrometer) spectrum (5 cm -1 resolution), obtained 
at Palestine, Texas, shown in Fig. 1 reveals marked spectral structure in the 
region of 9.6/a ozone band. This feature present in observed spectrum prompted 
the present investigation to indirectly determine the atmospheric ozone. 

The narrow (about 1.5/a wide) 9.6/a ozone band in such emission spectrum of 
the earth and atmosphere is influenced by the water vapor absorption in the so 
called 'window' region extending from 8 to 12/a. So it is necessary to know the 
distribution of water vapor in the atmosphere in order to interpret the spectrum 
in the S.6/a region. Although in principle it is possible to determine the distri- 
bution of water vapor, temperature and ozone indirectly from the radiation 
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measurements we have not attempted to do so in this preliminary study. Instead, 
the distribution of ozone alone is determined from the spectrum when the water 
vapor and temperature distributions are known. 

Earlier studies which attempted to determine the atmospheric ozone from 
ground based measurements of the emission spectrum of ozone in the 9.6/x region 
were not encouraging (see for e.g., Dave et al., 1963). However some disad- 
vantages encountered in interpreting the atmospheric emission spectra when 
measured at ground may not be encountered when the spectrum is obtained from 
space. 

With this view in mind we have computed spectra for several stations at 
which soundings of temperature, water vapor and ozone were available. These 
synthetic spectra were then subjected to an inversion procedure to deduce the 
ozone distribution. The inversion procedure adopted here depends on a two 
parameter model for describing the ozone distribution. Such a model repre- 
sentation is arrived at from broad based climatological features of the atmos- 
pheric ozone. 

The feasibility of determining atmospheric ozone is then demonstrated from 
an error analysis of four case studies. A small finite error in radiance is mag- 
nified to a large error in the inferred atmospheric ozone content. The magnifi- 
cation of error is different for each case and is intolerably large under certain 
circumstances. 
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On account of several inadequacies of data pertaining to the absorption 
coefficients and atmospheric temperature and water vapor sounding we could not 
perform an inversion and obtain the ozone distribution for Palestine, Texas. 
However a reasonable match of the spectrum is shown to illustrate the general 
soundness of the scheme. 

This investigation is intended to be a pilot study of this problem. Several 
spectra will be available from the forthcoming meteorological satellites which 
may be analysed with the help of the findings from this study. 

MODEL OF THE ATMOSPHERIC OZONE DISTRIBUTION 

The atmospheric ozone is observed to change appreciably from day to day 
below about 25 km. At higher elevations, where photochemical equilibrium is 
restored in a day or less, ozone changes are systematic and follow a seasonal 
variation (Dutsch, 1963). In Fig. 2 the mean meridional distribution of ozone for 
spring and fall are shown (after Hering and Borden, 1965, a). Some significant 
features, concerning the vertical ozone distribution, can be deduced from these 
meridional cross sections. The pronounced ozone maximum in the stratosphere 
is present at all latitudes. The position of ozone maximum which is at about 
24 km at low latitudes drops to lower altitudes as the north pole is approached. 
Also the rather sharp maximum present at low latitudes is replaced by a con- 
siderably broad one at high latitudes. The position of the tropopause as a func- 
tion of latitude has a strong influence on the meridional distribution of ozone. 

The strong stratification that is present in the stratosphere is replaced by a very 
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weak structure in the troposphere. On a further examination of several ozone 
soundings one can observe that the tropopause in the tropical latitudes in all 
seasons marks the level of transition from strong stratospheric to weak tropo- 
spheric structure in ozone. On the other hand at middle and high latitudes the 
tropopause has such significance only in summer and fall. During winter and 
spring such transition takes place at a level 3 to 4 km below the tropopause 
height. 

With these climatological features a model of the distribution of ozone is 
made involving minimum number of parameters to describe it. Photochemical 
equilibrium distribution is assumed above some high altitude Z x , (^30 km). 
Further down from the level Z to the tropopause height, Z 2 , a Gaussian dis- 
tribution, which can simulate the pronounced maximum is adopted. Then in the 
troposphere the ozone concentration is assumed constant having a value equal to 
that obtained at Z 2 , the tropopause height. 

Let the Gaussian distribution describing the ozone variation between Zj and 
Z 2 be given by 


0 3 (Z) 



(1) 


where 0 3 (Z) represents the ozone concentration at any height Z, (z x > Z > Z 2 ), 
and 0 3 (Z m ) is the ozone concentration maximum at the height z . h is a measure 
of the dispersion of ozone about the height Z . 
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Now applying Equation (1) at the selected high level Z t , where 0 3 (Zj) is 
assumed from Photochemical equilibrium, we have 


° 3 ( z i) " ° 3 ( z m) ex P ~ [' '' h"] 

From Equations (1) and (2) 0 3 (Z) may be eliminated to yield 


( 2 ) 


0 3 (Z) = 0, 
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(3) 


Equation (3) describes ozone distribution, in the height region Z l to Z 2 , in 
terms of two unknowns Z m and h . This is essentially the two parameter 
representation. 


THEORY 

We know from radiative transfer theory that intensity I , of radiation at 
frequency v, upwelling from the earth's atmosphere is given by the equation 
(see for e.g., Elsasser, 1960) 


I„ = B„ (T„K(0> + [ B v (T,)<tr„(Z) 

Jr „(°) 


(4) 


where 

B v is the Plank intensity, 

T 0 and T z are the absolute temperatures at the ground and any height Z, 
r v (0) and r v (Z) are the transmissions from ground and any level Z to top 
of the atmosphere. 
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The Equation (4), which forms the basis of all indirect methods or 'inver- 
sions,' is applied to small but finite, 5 cm -1 wide (to correspond to IRIS resolu- 
tion), spectral intervals of the emission spectrum at 9.6/u. region. As the 9.6p 
band of ozone is about 150 cm -1 wide we may think of having about 30 such use- 
ful spectral intervals yielding as many independent pieces of information. How- 
ever in practice we find that only two independent pieces of information can be 
elicited. This point will be elaborated in the discussion of the results. 

The transmission r v in the 9.6/x region is influenced by both water vapor 
and ozone. So the effective transmission is given by 

r v = ( t h 2 o ' r o 3 ) v (5) 

where t h 0 and r Q are the transmissions of water vapor and ozone. 

Let us choose two spectral intervals labled 1 and 2, in the 9.6/i band, for 
which the upwelling intensities may be written as 

I, = B, (T 0 )r 1 (0) ♦ Jb, (T,)dT, (Z) (6) 

= B 2 (T„)t 2 (0) + Jb 2 (T,)^ (Z) (7) 

Since we are assuming the vertical soundings of temperature and water 
vapor we need to solve for the two parameters describing the vertical ozone 
distribution from Equations (6) and (7) . The information pertaining to ozone is 
contained in the transmission functions v and r 2 . 
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Let us suppose a small increment SI in the intensity be due to a small in- 
crement SZ m in the position of ozone maximum and a small increment Sh in the 
dispersion. We may express such relationships by the equations 



<9I 2 <9 1 2 

SI 2 = dl Sz m + “dlT sh ( 9 ) 

m 

The above simultaneous equations may be used to correct, by an iteration 
procedure, some initial guess for Z m and h and eventually obtain correct values 
of these parameters, describing the ozone distribution, that yield the desired 
intensities and I 2 . This is in essence the perturbation procedure used here. 


TRANSMISSION FUNCTION 

Based upon several earlier studies Moller and Raschke (1963) have presented 
a method of obtaining the transmission function for the water vapor in the 8-12p, 
region. This method is adopted in the present study. 

The transmission function for ozone is taken from a study of Walshaw (1954), 
where he has investigated the 9.6p band with a resolution of 6.5 cm - *. Walshaw 
has fitted the transmissio' oy means of a statistical band model with allowance 
for a ’modulated contour’ for the band (see Goody, 1964, a). The transmission 
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function is given by 


where 


sinh k 1/2 X 


V V 


exp " (Xj] 


kj , 172 


sin k 


1/2 






when k > 0 


when k < 0 


> ( 10 ) 



< j v is the mean line strength* 
a is the mean line half width 
8 is the mean line spacing 
m is the path length of the absorbing gas 
k v is an emperically determined constant as a function of v to get the best 
fit to the experimental measurements. 


The above relationship is based on the assumption of exponential intensity dis- 
tribution of spectral lines. 

Walshaw's measurements were made for homogeneous paths in the laboratory 
while in the atmosphere we do not obtain such homogeneity conditions as the ab- 
sorbing gas is distributed along a strongly varying pressure path. The variation 


* All the parameters cr v , a, Sand are obtained from Walshaw (1954). However according to 
Walshaw (1957) a was multiplied by 1.30. 
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of temperature along the path though present in the earth's atmosphere is never 
as pronounced. To take into account such inhomogeneous conditions, the Curtis- 
Godson approximation is commonly used. This approximation is shown to be 
quite inadequate for the atmospheric ozone by Walshaw and Rodgers (1963) 
amongst others. It is possible to incorporate such inhomogeneous conditions 
and obtain radiative transfer without any approximation in a single spectral line 
(see for e.g., Hitschfeld and Houghton, 1961). Now the transmission r' in a 
spectral interval, large compared to a line, may be expressed in terms of the 
average absorption in one line provided we adopt a statistical model, for the 
transmission function, which assumes lines of equal intensity but with random 
position (Goody, 1964, b). Thus we have 

r' = exp - -g- { H) 

where W^, the equivalent width of a line, may be calculated numerically from 



Here 

a is Lorenz half width 
v the frequency at the center of the line. 

Plass (1938) has shown that one can match this statistical model of the 
transmission function having uniform line intensity with a statistical model 
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having an exponential intensity distribution for lines. From such a method of 
matching the two types of statistical band models we arrive at the following 
transmission function for the 9.6p ozone band. 
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exp 


- tv] 


sinh k 1/2 X ' 


k 1/2 X ' 

V V 


~\ 


when 


k > 0 
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(13) 


sin |kj 1/2 \ 


Ik. I 1/2 X. 


when 


k v < 0 
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where 


X 
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7T H 
4 TT 


and the other symbols have the same meaning as before. 


DISCUSSION OF RESULTS 

The main objective of this study is to examine the feasibility of determining 
atmospheric ozone. Since we do not at present have spectra taken at more than 
one station it has been necessary to calculate emission spectra from several 
known simultaneous soundings of ozone, temperature and dew point temperature, 
and then to invert such synthetic spectra to get back the ozone distributions. We 
have done this with the help of the atmospheric data presented in ’ozone sonde 
observations over North America’ by Hering and Borden (1965). Three such in- 
versions are shown in Figs. 3, 4 and 5. The temperature and ozone data were 
extrapolated to 50 km height. 
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Two spectral intervals, 5 cm -1 wide, centered at 1000 cm -1 and 1050 cm -1 
are used in the inversion scheme. These spectral intervals are chosen such that 
one of them at 1000 cm" 1 has moderate ozone absorption while the other at 
1050 cm" 1 has strong absorption. 

A comparison of the computed and observed distribution of ozone at Befora 
(Fig. 3) shows the ability of the computed curve to approximate smoothly the ob- 
served one. Although the computed curve can't reproduce the double maximum 
it succeeds in predicting the total ozone with 2% accuracy. This is so on ac- 
count of the extremely small convergence level, 0.1%, achieved in matching the 
radiance computed from observed and calculated ozone curves. In Fig. 4 we see 
the solution obtained for an ozone sounding at Fairbanks. Although the observed 
ozone profile in this case is not as complicated as the one shown for Bedford, 
the calculated total ozone is off by 5%. This is primarily because of the 1% con- 
vergence level to which the radiances were matched. The computed spectrum 
for each case is shown in the respective figures. The 9.6/x band in these two 
cases stands out quite significantly as there is a good amount of stratification 
in the temperature. To emphasize this aspect we show in Fig. 5 a case study 
for an Arctic Station, Thule, Greenland, where the temperature stratification is 
very weak resulting in almost a black body continuum for the spectrum. The 
inverted ozone profile, calculated with such a weak spectrum, shows a poor 
agreement with the observed profile. The computed total ozone is in error by 
50%. 
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The observed spectrum at Palestine, Texas, offers a means to test the in- 
version scheme. So an attempt was made to invert the observed radiances in 
the 9.6/;. region of the spectrum and obtain an ozone distribution. This attempt 
was not successful on account of two important reasons. The meteorological 
data needed for inversion, namely the temperature and dew point temperature 
soundings at Palestine were not available. So it was necessary to use such data 
from a neighboring meteorological station, Shreveport. Secondly the absorption 
coefficients adopted for water vapor and ozone, although reasonable, are not 
exact. It is very important to have accurate measurements of the absorption 
coefficients in order to have a successful inversion scheme. So as an alterna- 
tive to a rigorous inversion we have attempted to match the observed Palestine 
spectrum with the help of the Shreveport soundings and an ozone profile charac- 
teristic of the tropics shown in Fig. 6. The observed and computed spectra are 
compared in Fig. 7. 

At this point it is desirable to examine the question concerning the number 
of independent pieces of information about atmospheric ozone that one can possibly 
elicit out of the emission spectrum in the 9.6 p region. This may be demonstrated 
from the nature of the weighting functions. In Fig. 8 we show, for the case of 
Fairbanks, the weighting functions dr/dzfor ozone, in the two spectral regions, 
of 5 cm -1 width, centered at 1000 and 1050 cm" '. The weighting functions re- 
semble one another and essentially have the shape of the ozone profile. This is 
so on account of the weak absorption in the 9.6/j. region. Even in the interval at 
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1050 cm -1 , near the center of the absorption band, the ground is 'visible' and a 
significant part of the ground emission is transmitted to the top of the atmos- 
phere. From the strong similarity of the two weighting functions we may infer 
that they cannot convey substantially different information about the ozone dis- 
tribution. Nevertheless a small amount of dissimilarity that is present in the 
two weighting functions makes the information content of the two spectral intervals 
to some extent independent of one another. This is the reason for limiting our- 
selves to two independent parameters. An attempt to get three pieces of infor- 
mation with three spectral intervals resulted in a failure. 

Let us now examine the practicability of determining the atmospheric ozone 
indirectly from measured radiance in the 9 ,6p region. As there are in practice 
some errors present in the measured radiance, we have to examine the degree 
to which the deduced solution will be influenced by such errors. For this purpose 
we have made an analysis of the error introduced in the total ozone due to a 
known error in radiance. 

Let the total ozone be denoted by [o 3 ] and a small error in it by 5 [o 3 ]. 
Similarly let I and 51 be the corresponding intensity and an error in it. Now 
since both the total ozone and intensity are functions of the two parameters Z m 
and h we seek to get 

/AlogJoJ / A log I \ 

V A l0gZ m / A l °e Z J. 
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which gives an amplification of the error in total ozone due to an error in in- 
tensity when h is held constant. In an analogous way we can get 

/ A log [o 3 ] /a jpg A 

\ Alogh j A log hy z 

the amplification of error when Z m is a constant. With the help of the four cases 
presented in Figs. 3 to 6 we have obtained these amplifications and they are 
tabulated for each case and the two wavelengths in Table 1. We find from this 
study there is generally an amplification of the error. An error of 1% in the 
measured radiance may result in an error of 5% in the derived total ozone. In 
the particular case of Thule, where the atmosphere is close to isothermal con- 
dition, the amplification of error is enormous and no meaningful information on 
ozone could be extracted. 

The above analysis did not consider the errors introduced by the assump- 
tions made in arriving at a model for the ozone distribution. The assumption of 
photochemical equilibrium above some high altitude, 30 km, and the height of 
the level separating the stratospheric and trophospheric ozone should be ex- 
amined. From Fig. 8 where the weighting functions are shown we see that as 
much as 20% of atmospheric emission could come from the region above 30 km 
and about 10% from the troposphere. By assuming photochemical equilibrium 
above 35 km we may be able to reduce the emission of the top layer by half. 
Even then small errors arising from these assumptions will be magnified in 
the solution obtained, as shown earlier. 
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There can be errors introduced from the incorrect temperature and dew 
point temperature soundings. Since about 50% of the emission from the ground 
is transmitted to the top we need to know the ground temperature and emissivity 
quite accurately. 

From the foregoing study we can arrive at the following conclusions. At the 
stations where radiosonde data are available total ozone could be estimated to 
about 10 to 20% accuracy provided the satellite measured radiance is accurate 
to 1% in the 9.6/x region. Although the accuracy of the indirect determination is 
poor we may hope to get at least a crude picture of the day to day variation of 
total ozone over the mid latitudes where 10 to 20% changes in total ozone occur. 
The indirect ozone determination may be improved by improving the spectral 
resolution of the radiance measured. 
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Table 1 


Amplification of Error 


Case 


d log [° 3 ] d lo e z m 

d log [0 3 ] $ logh 


d log Z m d log I 

d log h d log I 


I, 

-3.1 

-4.9 

FAIRBANKS 


*2 

-2.0 

-5.2 


I, 

-3.7 

-4.0 

BEDFORD 


l 2 

-2.6 

-3.3 


I, 

-34.5 

-105.0 

THULE 


*2 

-158.0 

+104.6 

SIMULATION 


-4.3 

+2.96 

OF PALESTINE 

i 2 

-1.5 

+0.4 
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igure 1. Emission Spectrum Measured at Palestine, Texas, May 8, 1966 
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LEGENDS TO FIGURES 

Fig. 1. Emission spectrum measured at Palestine, Texas, May 8, 1966. 

Fig. 2.a. Mean ozone density (u gms/m 3 ) for March, April, and May (1963, 

1964). After Hering & Borden 1965. 

b. Mean ozone density (u gms/m 3 ) for September, October, and November 
(1963, 1964). After Hering & Borden 1965. 

Fig. 3. Case Study: Bedford, Feb. 20, 1964. 

Fig. 4. Case Study: Fairbanks, June 24, 1964. 

Fig. 5. Case Study: Thule, March 25, 1964. 

Fig. 6. Emission spectrum calculated to match the spectrum obtained at 
Palestine on May 8, 1966. 

Fig. 7. Comparison of the Palestine emission spectrum in the 9.6/x region 
with the calculated spectrum shown in Fig. 6. 

Fig. 8. Weighting functions for ozone calculated for the case of Fairbanks, 
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BEDFORD FEB 20, 64 

TOTAL 0 3 = 0.362 Cm S.T.P. 
LEVEL OF CONVERGENCE 0.1% 
ERROR IN COMPUTED TOTAL 0 3 2% 



Figure 3. Case Study: Bedford, Feb. 20, 1964 





THULE MAR 25, 64 

TOTAL 0 3 = 0.550 cm S.T.P. 

LEVEL OF CONVERGENCE 1% 
ERROR IN COMPUTED TOTAL 0 3 + 50% 


E 



OZONE CONCENTRATION (cm S.T.P./Km) TEMPER. 



WAVE NO cm" 



26 


Figure 7. Comparison of the Palestine Emission Spectrum in the 9.6/x Region With the Calculated Spectrum Shown i 



WEIGHTING FUNCTIONS FOR FAIRBANKS JUNE 24, 1964 
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Figure 8. Weighting Functions for Ozone Calculated for the Case of Fairbanks, Alaska 



